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1. Introduction 

Quantum Chromodynamics(QCD) is believed to be the fundamental theory of strong inter- 
action. At finite temperature, QCD is usually described by two extreme pictures. One is with the 
weakly interacting meson gas in the low temperature regime and another is with perturbative Quark 
Gluon Plasma (QGP) in the high temperature regime. Lattice studies of QCD Equation of State 
(EOS) studies [|l|] indicate that the Steven-Boltzman ideal gas limit can be reached only at high T, 
and deconfined partons might be strongly interacting in the intermediate temperature range above 
T c . This point is supported both by RHIC experiments and theoretically studies. On the one hand, 
the QGP observed by RHIC can be well described by the hydrodynamical model [Q]. On the other 
hand, the lattice studies of meson correlators show that charmonia and light hadrons can survive in 
a temperature range beyond T c [Q, [|, Bp. 

Since glueballs are predicted by QCD and well defined in pure Yang-Mills theory, their T- 
evolution is also a good probe to investigate the property of QCD matter in the deconnnement 
phase. In this work, we carry out a numerical lattice study on anisotropic lattices with much 
finer lattice in the temporal direction than in spatial ones. By varying the temporal extension of 
the lattice, we obtain a wide temperature range from 0.3T C to \.9T C . At each temperature, we 
take into account all the twenty R PC channels, with PC = ++,H — , — h, the various parity- 
charge conjugate and/? =A\,Az,E,Tt,T2 the irreducible representations of lattice symmetry group. 
For each R PC , as in the zero temperature case [Q, [7|, ||], we implement smearing schemes and the 
variational method to acquire an optimal glueball operator which couples most to the ground state. 
In the data analysis stage, the correlators of these optimized operators are analyzed through two 
approaches. First, the thermal masses Mq of glueballs are extracted in all the channels and all over 
the temperature range by fitting the correlators with a single-cosh function form, as is done in the 
standard hadron mass measurements. Thus the T-evolution of the thermal glueball spectrums are 
obtained. Secondly, with the respect that the finite temperature effects may result in mass shifts and 
thermal widths of glueballs, we also analyze the correlators in A\ + , A\ + , E ++ , and T 2 ++ channels 
with the Breit-Wigner ansatz which assumes these glueballs thermal widths, say, changes Mq into 
a\) — iT in the spectral function (see below). It is expected that the temperature dependence of (Oq 
and r can shed some light on the scenario of the QCD transition. 

2. Numerical details 

We adopt the tadpole improved Symanzikars action, which has been extensively used in the 
study of glueballs, 



where j8 is related to the bare QCD coupling constant, t, = a s /a t is the aspect ratio for anisotropy 
(we take E, = 5 in this work), u s and u t are the tadpole improvement parameters of spatial and 
temporal gauge links, respectively. Lattices are 24 3 x N t with j8 = 3.2, and a s = 0.0878fm, and the 
spatial volume V ~ (2.1fm) 3 . 
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The deconfinement critical point is roughly determined by the susceptibility %p of Polyakov 
loops, 

Xp = (® 2 ) - (@) 2 (2.2) 
where denotes the Z(3) rotated Polyakov line, 

{RePexp[-27ri/3]; argP e [7f/3,7f) 
ReP; argPG [-71/3,71/3) . (2.3) 

RePexp[27r//3]; argP e [-n, -n/3) 

As illustrated in Fig. |IJ, the peak position is roughly at N t = 38, which corresponds to the 
critical temperature T c . With the lattice spacing a s = 0.0878fm, T c is estimated to be T c = 296MeV. 
Based on this, we vary N t to obtain a temperature range from T ~ 0.3P ( to 1 .9T C , as listed in Table |]. 
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Figure 1: %p * s plotted versus at /3 = 3.2. There is a peak of %p near N, = 40. 

For all the 20 J PC channels of glueballs, we take the following two steps to construct the 
optimal glueball operators which couple most to the ground states (More details can be found in 
Ref. [Q] and Ref. []7j, |8|]). First, for a given gauge configuration, we generate six differently smeared 
copies, on each of which, four realizations of each J PC are established based on all the different 
spatially oriented Wilson loops of a set of loop prototypes [0]. As a result, we obtain a set of 24 
different glueball operators, = {(j) a , a = 1,2, .. . ,24}, for each J PC . At the second step, we carry 
out the variational method on each operator set (p to determine the specific combinational coeffi- 
cients {v a ,cc = 1,2, .. . ,24} relevant to the ground state, such that the desired optimal operator is 
obtained as <1> = £v a a . Practically at each temperature, after a thermalization of 10000 heatbath 
sweeps, the glueball operators are measured every three compound sweeps, each of which is com- 
posed of one heatbath and five micro-canonical over-relaxation(OR) sweeps. In the data analysis, 
the measurements are divided into A^bin bins of the size n m \, = 400 . Parameters A^m and n m b at 
various temperature are also listed in Table [j] 

3. Glueballs at finite temperature 

Theoretically, under the periodic boundary condition in the temporal direction, the temporal 
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Table 1: Temperature range and simulation parameters in this work. 



correlators C(t,T) at the temperature T can be written in the spectral representation as 
C{tJ)^^-Tv(e- H l T <S>{t)<S>{0)) 



Z(T) 

CO 

= J dcop{co)K(co,T), 



' 2T ^ ^ n ^" 



(3.1) 



with a T-dependent kernel K(co, T) = ^^{mjJiT))^ an( * tne s P ectra l function, 

P(<0) = I Kn ^l™l l2 e- EjT (S(co - (E n - E m ) - 8(co - (E m - £„)), (3.2) 

m,n ) 

where Z{T) is the partition function at T, and E n the energy of the thermal state \n) (|0) represents 
the vacuum state). In the zero-temperature limit(r — > 0), due to the factor exp(—E m /T), the spectral 
function p (w) degenerates to 



p(«)=E 



2Z(0) 



(8((0-E n )-8((0 + E n )), 



(3.3) 



thus we have the function form of the correlation function, C(t,T = 0) = £W„e E "' with W n = 

n 

|(0|<!>|«}| 2 /Z(0). However, for any finite temperature (this is always the case for finite lattices), all 
the thermal states with the non-zero matrix elements (m\<t>\n) may contribute to the spectral func- 
tion p{(o). Intuitively in the confinement phase, the fundamental degrees of freedom are hadron- 
like modes, thus the thermal states should be multi-hadron states. If they interact weakly with each 
other, we can treat them as free particles at the lowest order approximation and consider E m as the 
sum of the energies of hadrons including in the thermal state |m). Since the contribution of a ther- 
mal state \m) to the spectral function is weighted by the factor e\p(—E m /T), apart from the vacuum 
state, the maximal value of this factor is exp(— M m ; n /r) with M m [ n the mass of the lightest hadron 
mode in the system. As far as the quenched glueball system is concerned, the lightest glueball 
is the scalar, whose mass at the low temperature is roughly M ++ ~ 1.6 GeV, which gives a very 
tiny weight factor exp(— M Q ++/T C ) ~ 0.003 at T c in comparison with unity factor of the vacuum 
state. That is to say, for the quenched glueballs, up to the critical temperature T c , the contribution 
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of higher spectral components beyond the vacuum to the spectral function are much smaller than 
the statistical errors (the relative statistical errors of the thermal glueball correlators are always a 



few percents) and can be neglected. As a result, the function form of p (ft)) in Eq. [O] can be a good 
approximation for the spectral function of glueballs at least up to T c . Accordingly, considering the 
finite extension of the lattice in the temporal direction, the function form of the thermal correlators 
can be approximated as 

cosh(M„(l/(2r)-f)) 



(3.4) 



sinh(M„/(2r)) ' 

which is surely the commonly used function form for the study of hadron masses at low tempera- 
tures on the lattice. As is always done, the glueball masses M n derived by this function are called 
the pole masses in this work. 

3.1 The single-cosh fit 

After the thermal correlators C(f, T) are obtained, the pole masses of the ground state (or the 
lowest spectral component) can be extracted straightforwardly. For each R PC channel and at each 
temperature T, we first calculate the effective mass M e ff (t) as a function of t by solving the equation 



C(t + \,T) cosh((f + 1 -N t /2)a t M en {t)) 



C(t,T) 



cosh((7 -N t /2)a t M e ff(t)) 



(3.5) 



and determine the time window [?i,?2] where M e ff(t) has a plateau. In this time window, C(t,T) is 
fitted through a single-cosh function form. Fig. illustrates this procedure in E ++ channel. 
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Figure 2: Effective masses at different temperatures in E ++ channel. Data points are the effective masses 
with jackknife error bars. The vertical lines indicate the time window [fi,?2] over which the single-cosh 
fittings are carried out, while the horizonal lines illustrate the best fit result of pole masses (in each panel the 
double horizonal lines represent the error band estimated by jackknife analysis) 



In this work, the pole masses in all the 20 R PC channels are extracted at all the temperatures. 
The common feature of temperature dependence of pole masses is that, below T c , the pole masses 
keep stable when varying the temperature, while above T c , the pole masses decrease gradually and 
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Figure 3: The T-dependence of pole masses Aj , A^ + and T 2 ++ glueballs. 



cannot be extracted beyond the temperature T = \.6T C . Fig. ^illustrates these hebaviors in A\ + , 
A^ + , E ++ , and T 9 ++ channels. These results are different from the observation of a previous lattice 



study on glueballs where the observed pole-mass reduction start even at T ~ 0.8r f [10]. 



3.2 Breit-Wigner analysis (BW) 

Theoretically in the deconfined phase, gluons can be liberated from hadrons. However, the 
study of the equation of state shows that the state of the matter right above T c is far from a pertur- 
bative gluon gas. In other words, the gluons in the intermediate temperature above T c may interact 
strongly with each other and glueball-like resonances can be possibly formed. Thus different from 
bound states at low temperature, thermal glueballs can acquire thermal width due to the thermal 
scattering between strongly interacting gluons and the magnitudes of the thermal widths can signal 
the strength of these type of interaction at different temperature. 



By assuming glueballs thermal widths, we also adopt the Breit-Wigner ansatz [jlOJ to analyze 
the thermal correlators once more. First, we treat thermal glueballs as resonance objects which 
correspond to the poles (denoted by ft) = (Oq — iT) of the retarded and advanced Green functions in 
the complex ft)— plane, coq is called the mass of the resonance glueball and T its thermal width in 
this work. Secondly, we assume that the spectral function p(ffl) is dominated by these resonance 
glueballs. With these respects, the thermal correlators can be parameterized as [R [T0|] 



grit) 



Re 



cosh((ojo + /r)(2^ 



0) 



sinh 



(<BD+'T) 



2T 



+ 2(OqT £ cos (InnTt) 



n=\ 



(2nnT + T) 2 + w 2 



(n 



(3.6) 



In the data analysis, the effective mass ft) ( , (?) and effective width T (?) are obtained first by 
solving the equation array, 



gr(t)/gr(t + l) =C(t,T)/C(t + l,T) 
gr(t + l)/gr(t + 2)= C(? + 1, T)/C{t + 2, T) 



(3.7) 
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Figure 4: Determinations of the fit range [fi,?2] in r 2 ++ channel at N t = 128, 36, and 20. In each row, 
®0 (0 anc ^ r^ eff '(f) are plotted by data points with jackknife error bars. [t\,t2] are chosen to include the 
time slices between the two vertical lines, where COq^ (f ) and r' eff ) (/) show up plateaus simultaneously. The 
best fit results of coq and F through the function gr(t) are illustrated by the horizonal lines. 



where C(t,T) is the measured correlator, then the simultaneous plateau region of (Oo(t) and T(t) 
gives the fit window \t\ , t<z\ where the fit is carried out. This procedure in T 2 ++ channel is shown in 
Fig. [4] for instance. 

The main feature of the best fit (Oq and T in Af + ,A^ + ,E ++ , and T 2 ++ channels is illustrated 
in Figure ||: First, (Hq's are insensitive to T, or more specifically, the reduction of &>o at the highest 
temperature T = 1.90r c are less than 5%. Secondly, T's are small and do not vary much below T c , 
but increase abruptly when the temperature passes T c and reach to values ~ (Oq/2 at T = 1.90r c . 
These features can be easily seen in Fig. where the behaviors of o\) and T with respect to the 
temperature T are plotted for all the four channels. 
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Figure 5: £0o's and T's are plotted versus T /T c for Aj~ + , Aj + and T 2 ++ channels. The vertical lines 
indicate the critical temperature. 



4. Summary and conclusion 

In the pure SU(3) gauge theory, the thermal correlators in all the 20 symmetry channels are 
calculated on anisotropic lattices in the temperature range from 0.3 T c to 1.9r f . Both the single- 
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cosh fit and BW analysis show that glueballs can survive up to 1.97^. Our results are consistent 
with that of the studies of EOS and charmonia |], 12, [l3|]. In BW analysis, glueball masses 
keep stable when the temperature increasing, and the thermal widths of glueballs becomes larger 
and larger above T c . It seems that in the intermediate T range, the state of matter are dominated by 
strongly interacting gluons. Gluons interact with each other strongly enough to form glueball-like 
resonances, in the mean time, glueballs can also decay into gluons. At a given temperature, these 
two procedure reach the thermal equilibrium. The thermal widths signal the interaction strength. 
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